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ABSTRACT 


Results  of  a  computer  study  which  theoretically  determined  the 
detection  and  location  capabilities  of  various  wavenumber  spectra  techniques 
are  presented  in  thic  report.  Two  distinct  short-period  array  configurations 
were  simulated.  The  first  configuration  simulated  was  a  13- element  array 
employing  the  outputs  of  the  inner  13  subarrays  of  the  Montana  LASA.  In 
this  case,  the  ambient  noise  was  considered  to  be  completely  uncorrelated. 

The  second  configuration  simulated  consisted  of  all  or  a  portion  of  a  standard 
LASA  subarray.  Various  noise  models,  made  up  of  mostly  organized  noise, 
were  used.  Four  spectral  techniques  —  the  conventional  ■  pectrum,  two 
high-resolution  spectra,  and  the  probabilistic  processor  —  were  evaluated. 
Unfortunately,  an  error  in  the  program,  which  could  not  be  rectified  in  time 
for  inclusion  of  correct  results  in  this  report,  prevented  evaluating  the  de¬ 
tection  capability  of  the  probabilistic  processor. 

When  wavenumber  spectra  at  a  single  frequency  are  used  to 
locate  a  signal  in  uncorrelated  noise,  all  four  techniques  must  yield  exactly 
the  same  location.  It  was  found  that  events  with  an  SNR  of  1.  0  or  greater 
were  almost  always  locaced  properly  within  the  resolving  power  of  the  grid 
points  at  which  the  spectra  were  calculated.  Wideband  frequency-domain 
location  can  be  accomplished  by  computing  wavenumber  spectra  at  several 
independent  frequencies  and  summing  the  spectra  in  a  velocity-preserving 
manner  before  choosing  the  peak-power  location. 

Conventional  spectra  of  signals  in  uncorrelated  noise  were 
computed  at  five  independent  frequencies  and  summed.  This  procedure  re¬ 
sulted  in  no  incorrect  locations  for  SNR's  of  0.  5  or  greater,  thereby  indicating 
its  significant  superiority  over  the  location  ability  of  single-frequency  spectra. 
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The  wideband  wavenumber  spectral  technique  is  roughly  analogous  to  time- 
domain  beamforming,  and  the  results  obtained  here  may  be  indicative  of 
those  achievable  in  the  time  domain. 

The  performance  of  the  small  arrays  in  organized  noise  was 
found  to  be  markedly  inferior  to  the  performance  of  the  large  array  using 
uncorrelated  noise.  The  conventional  spectrum  seems  to  have  no  practical 
ietection  ability,  while  the  high-resolution  techniques  are  perhaps  marginal 
in  this  respect.  All  techniques  located  a  significant  number  of  events  im¬ 
properly,  even  when  the  SNR  was  as  high  as  6.  48. 
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SECTION  I 


INTRODUCTION 

The  relationship  oetween  the  detection  or  location  of  short- 
period  seismic  events  using  wavenumber  spectra  and  the  signal-to-noise 
ratio  (SNR)  of  the  data  being  processed  was  studied,  and  the  results  are  pre¬ 
sented  in  this  report. 

The  need  for  an  effective  detection  mechanism  to  monitor 
underground  nuclear  tests  is  obvious.  Reasonably  accurate  location  techniques 
assist  in  the  classification  problem,  since  events  from  many  areas  can  be 
rejected  immediately.  More  importantly,  the  approximate  location  of  an 
event  is  required  for  effective  off-line  beamsteer  or  MCF  processing  per¬ 
formed  to  enhance  classification  parameters. 

This  study  was  initially  undertaken  analytically  but  soon  became 
intractable.  The  problem  was  then  divided  into  two  experiments  designed  to 
determine  empirically  the  location  and  detection  abilities  of  particular  array 
geometries  as  a  function  of  SNR.  The  two  tests  treat  correlated  noise  and  un¬ 
correlated  noise  separately  and  process  correlated  noise  at  the  subarray  level 
and  uncorrelated  noise  using  LASA  A-,  B-,  C-,  and  D-ring  subarrays.  This 
separation  into  two  experiments  results  from  the  increased  computational 
difficulty  in  contrasting  correlated  noise  to  uncorrelated  noise.  Both  data 
generation  and  processing  techniques  require  matrix  rather  than  vector 
arithmetic  when  the  noise  field  is  correlated. 

For  each  of  the  two  experiments,  four  techniques  were  evalu  ¬ 
ated,  with  all  processing  in  the  frequency-wavenumber  domain.  The  tech- 

I 

niques  are  the  conventional  wavenumber  spectrum,  the  probabilistic  processor, 
and  two  forms  of  high-resolution  wavenumber  spectra.  Both  forms  of  high- 
resolution  wavenumber  spectra  were  found  to  have  equivalent  averaged  prop- 
erties  for  a  uniform  plane -wave  signal  and  normalized  f-k  spectra;  there¬ 
fore,  this  report  usually  reflects  the  averaged  properties  of  both  techniques. 
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For  the  processing,  a  random-number  generator  generates  a 
sequence  of  normal  (mean  =  0,  variance  =  1)  random  numbers  which  are 
paired  and  grouped  to  form  a  set  of  100  random  complex  vectors.  The  set 
of  vectors  is  transformed  to  a  second  set  of  random  complex  vectors  having 
a  specified  covariance  matrix.  The  latter  vectors  are  taken  to  be  noise 
transform  vectors  with  the  specified  noise  crosspower  matrix.  A  scaled 
signal  transform  vector  is  added  to  each  noise  transform  vector  to  form 
the  data  transform  vector.  The  signal  is  scaled  to  yield  a  desired  signal- 
to-noise  power  ratio.  Resulting  data  are  then  processed  by  each  of  the  pro¬ 
cessors  at  selected  points  in  the  vicinity  of  the  signal  location  in  K-space. 

In  the  correlated  noise  cases,  the  values  obtained  for  each  noise  vector 
over  K-space  are  normalized  to  have  a  unit  sum.  The  mean  and  standard 
deviation  at  each  wavenumber  point  for  each  processor  are  the  outputs  of 
the  computer  program.  The  program  steps  may  be  repeated  using  the  same 
signal  and  noise  vectors  but  with  a  different  signal  scale  factor.  An  error 
was  detected  in  the  program  which  generates  this  output  for  the  probabilistic 
processor.  Since  time  did  not  permit  correction  of  the  error,  the  means 
and  standard  deviations  are  shown  only  for  the  conventional  and  highrresolution 
schemes. 

The  above  steps  yield  mean  wavenumber  spectra  for  each  of 
the  four  processing  techniques  at  various  signal- to-noise  ratios  and  various 
noise  models.  From  these  results,  it  is  possible  to  observe  the  mean  be¬ 
havior  of  a  detection  scheme  as  a  function  of  SNR  for  the  various  types  of 
noise  fields. 

A  second  type  of  output  was  also  generated.  For  each  of  the 
100  signal-plus-noise  vectors,  the  location  of  peak  power  in  K-space  is 
determined.  These  results  are  then  plotted,  showing  the  number  of  times  out 
of  100  that  the  peak  power  occurs  at  various  K-space  locations.  This  type  of 
display  provides  an  indication  of  the  confidence  that  can  be  attached  to  the  peak- 
power  location  for  various  SNR's. 


-2 


•donee  services  division 


SECTION  II 


PROCEDURES 


A.  PROCESSING  EQUATIONS 

The  four  schemes  and  the  corresponding  processing  quantities 

evaluated  are 

•  Conventional  Wavenumber  Spectrum  (C3) 


V 


*T 


V 


Probabilistic  Processor  (PP) 
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|n  +  vv*Tl 


exp 


*T 

N  +  VV 


-1 


j  x*T  In  +  vv  *1  x 


E — 1 
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High-Resolution  Wavenumber  Spectrum  (HR1) 


v*Tv - v*Txx*Tv 


6+  X  X 


•  High-Resolution  Wavenumber  Spectrum  (HR2) 

2 


.  *T 

«T  26  +  X  X_  y*Txx*Tv 
(6  +  X*TX)2 


In  these  equations, 

N  is  the  noise  covariance  ma  rix 
X  is  the  vector  of  data  transforms 
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V  =  V(k)  is  the  vector  of  phase  shifts  appropriate  for  the  frequency 
and  velocity  of  a  plane  wave  propagating  across  the  array 


|  |  denotes  the  determinant  of  the  enclosed  matrix 
6  is  a  scalar 


*  is  conjugate 
T  is  transpose 

Brief  derivations  of  the  probabilistic  processor  and  high- 

resolution  equations  are  given  in  Appendix  A.  Appendix  B  provides  an 

i  *Ti 

efficient  method  for  evaluation  of  |N  +  VV  |. 

For  the  case  in  which  the  noise  is  uncorrelated,  the  covariance 
matrix  N  is  the  identity  matrix  I,  and  the  data  manipulations  reduce  to  vector 
arithmetic  for  both  data  generation  and  data  processing.  Appendix  B  shows 

that 

|n  +  vv  *T  |  =  In  |(i  +  s*VTvs) 


where 

*T 

SNS  =  I 

I  =  identity  matrix 

For  N  si,  the  determinant  simply  becomes  1  +  M,  where  M  is  the  number 
of  channels. 
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*T 


The  quadratic  form  in  the  PP  equation  may  be  written  as 
-1 
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For  N  =  I,  this  reduces  to 


*t  *T 

*T  X  VV  X 
X  X  -  * 


*T 

1  +  V  V 


The  term  X*TX  is  the  sum  of  the  power  in  each  channel  at  the 
frequency  being  processed,  and  V  V  =  M,  the  number  of  channels.  Thus,  using 
uncorrelated  noise,  the  quadratic  forms  of  the  high-resolution  schemes  and 
the  probabilistic  processor  scheme  are  different  combinations  of  the  con¬ 
ventional  wavenumber  spectrum,  the  sum  of  the  power  in  each  channel  and 
the  number  of  channels.  All  four  techniques  have  their  peak  values  at  the 
same  wavenumber  location.  Using  correlated  noise,  the  probabilistic  pro¬ 
cessor  does  not  necessarily  have  its  peak  at  the  same  location  as  the  other 
three  schemes. 

B.  DATA  GENERATION 

Let  N  denote  a  noise  crosspower  matrix.  It  is  desired  to 

generate  a  sequence  of  random  complex  vectors  ^i2 . ^im) 

having  the  covariance  matrix  N.  Let  denote  a  set  of  complex  vectors 

satisfying 

EY.  =  0 

i 

*T 

EY.  Y.  =  I 

i  i 

where  E  is  the  expected  value.  There  exists  a  matrix  S  such  that 

*T 

N  =  SS 


Then  Z  is  generated  by 


Z.  =  . 

i  i 
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Note  that 


EZ.  =  ESY. 
i  1 

=  SEY. 


=  0 


and 


*T 

EZ.Z. 

i  i 


sfeT  *T 
ESY.Y.  S 

l  l 


*T  *T 
=  SEY.  Y.  S 

i  i 


=  SIS 
=  N 


*T 


Let  s  be  a  unit  plane-wave  signal  transform  vector  and  form 

the  vector  X.  by 
1  1 

X.  =  Z.  +  as 

i  i 

where  a  is  a  scale  factor.  The  vector  X.  is  then  the  data  transform  vector. 
C.  UNCORRELATED  NOISE  PROCESSING 


In  the  experiment  using  uncorrelated  noise,  each  processor 
used  data  having  an  infinite  velocity  signal  from  13  channels  of  the  LASA 
A-,  B-,  C-,  and  D-ring  subarray  locations.  Figures  II- 1  and  II -2  show  the 
mean  and  standard  deviation  of  the  processor  outputs  as  a  function  of  wave- 
number  for  SNR's  of  0.  72  and  1.  12.  These  statistics  were  estimated  at  11 
equally  spaced  points  along  the  N-S  axis  from  the  origin  to  k  (wavenumber)  = 
0.  1  cycle/ km.  Figure  II -3  gives  the  mean  and  standard  deviation  as  a 
function  of  SNR. 
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In  addition  to  calculating  the  mean  and  standard  deviation, 
each  of  the  100  signal-plus -noise  vectors  was  processed  to  determine  the 
peak-power  location.  As  mentioned  previously,  each  of  the  four  techniques 
leads  to  the  same  location. 

Figure  II-4  gives  the  relative  frequency  of  occurrence  of  peak 
power  at  25  wavenumber  locations  equispaced  along  the  N-S  axis  between  the 
origin  and  0.  096  cycle/km.  These  computations  are  performed  at  a  single 
frequency  for  three  different  values  of  SNR.  Smoothing  can  be  achieved  by 
computing  wavenumber  spectra  at  several  independent  frequencies  within  the 
range  where  the  peak  SNR  occurs.  A  velocity-preserving  stack  of  the  indi¬ 
vidual  spectra  is  then  formed  before  choosing  the  peak-power  location.  This 
procedure  with  conventional  spectra  is  analagous  to  forming  time -domain 
beams  for  the  various  wavenumber  locations  and  choosing  the  largest  beam- 
steer  output. 

Figure  H-5  shows  the  relative  frequency  of  occurrence  after 
stacking  over  five  frequencies.  The  time  gate  used  in  obtaining  transforms 
from  real  data  must  be  short  in  order  to  contain  only  that  portion  of  the 
record  where  the  signal  is  largest  but  must  be  long  enough  to  provide  reason¬ 
able  crosspower  estimates  for  the  widely  spaced  seismometers.  For  short- 
period  processing  out  to  the  LASA  D-ring,  about  15  sec  is  a  good  compromise. 
From  sampling  theory,  this  results  in  an  independent  frequencv  increment 
of  1/15  =  0.  067  Hz.  Thus,  by  smoothing  wavenumber  spectra  over  five  adja¬ 
cent  independent  frequencies,  an  effective  band  of  0.  33  Hz  is  covered. 
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D.  CORRELATED  NOISE  PROCESSING 


Three  general  types  of  noise  fields  were  generated: 

•  An  isotropic  disk  centered  at  the  origin 
with  a  radius  of  0.  1  cycle/km  and  having 
20  percent  uncorrelated  noise  added. 

•  An  isotropic  disk  centered  at  k  =  0.  05 
(on  N-S  axis)  with  a  radius  of  0.  01  cycle/ 
km  and  having  1-percent  uncorrelated 
noise  added 

•  Iatropic  disks  centered  at  k  =  0.  05  (on 
N-S  axis)  with  a  radius  of  0.  025  cycle/ 
km  and  having  1-percent,  10-percent, 
and  20-percent  uncorrelated  noise  added 

The  first  two  models  used  a  standard  LASA  subarray  and  had 
25  channels  of  data.  The  last  model  used  the  center  and  outer  15  locations 
of  a  standard  LASA  subarray  and  had  16  channels  of  data.  All  processing  used 
an  infinite -velocity  signal.  The  equations  were  evaluated  at  11  equally  spaced 
points  on  the  N-S  axis  from  the  origin  to  k  =  0.  1  cycle/km.  Signal-to-noise 
ratios  ranged  from  about  0.  0004  to  0.  &.  Again,  the  three  types  of  displays 
shown  in  Figures  II-6  through  11-20  are  mean  and  standard  deviation  vs  wave- 
number,  mean  and  standard  deviation  vs  SNR,  and  relative  frequency  of 
occurrence  of  peak  power. 
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STANDARD  DEVIATION 


MEAN 


I 

HIGH  RESOLUTION,  MEAN  ^ 

HIGH  RESOLUTION,  STD  DEV 

CONVENTIONAL  SPECTRUM,  MEAN  I 


CONVENTIONAL  SPECTRUM,  STD  DEV 


WAVENUMBER  (cycle/km) 


Figure  II-2.  Mean  and  Standard  Deviation  Vs  k,  Uncorrelated  Noise,  SNR 


II -8 


science  services  division 


HIGH  RESOLUTION,  MEAN 
HIGH  RESOLUTION,  STD  DEV 
CONVENTIONAL  SPECTRUM,  MEAN 
CONVENTIONAL  SPECTRUM,  STD  DEV 


SIGNAL- TO -NOISE -RATIO 


Figure  II-3.  Mean  and  Standard  Deviation  Vs  SNR  at  k  =  0.  0  cycle/km, 
Uncorrelated  Noise 
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Figure  II-6.  Mean  and  Standard  Deviation  Vs  k,  Isotropic  Bodywave  Noise 
with  20-Percent  Uncorrelated  Noise  Added,  SNR  =  0.  0 
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Figure  II -8.  Mean  and  Standard  Deviation  Vs  SNR  at  k  =  0.  0  cycle /km,  Iso¬ 
tropic  Bodywave  Noise  with  20-Percent  Uncorrelated  Noise  Added 
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Figure  II -9.  Mean  and  Standard  Deviation  Vs  k,  Concentrated  Directional 
Noise  with  I -Percent  Uncorrelated  Noise  Added,  SNR  =  0.  60 


5 


science  services  division 


STANDARD  DEVIATION 


HIGH  RESOLUTION,  MEAN 
HIGH  RESOLUTION,  STD  DEV 
CONVENTIONAL  SPECTRUM,  MEAN 
CONVENTIONAL  SPECTRUM,  STD  DEV 


SIGNAL -TO- NOISE  RATIO 


Figure  II- 10.  Mean  and  Standard  Deviation  Vs  SNR  at  k  =  0.  0  cycle/km,  Con 
centrated  Directional  Noise  with  1-Percent  Uncorrelated  Noise 
Added 
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Figure  II- 11.  Mean  and  Standard  Deviation  Vs  SNR  at  k  =  0.  05  cycle/km,  Con 
centrated  Directional  Noise  with  1-Percent  Uncorrelated  Noise 
Added 
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Figure  11-12.  Mean  and  Standard  Deviation  Vs  k,  Directional  Noise  with  1 
Percent  Uncorrelated  Noise  Added,  SNR  =  0.  60 
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Figure  11-13.  Mean  and  Standard  Deviation  Vs  SNR  at  k  =  0.  0  cycle/km,  Direc 
tional  Noise  with  1-  Percent  Uncorrelated  Noise  Added 
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Figure  11-14.  Mean  and  Standard  Deviation  Vs  k,  Directional  Noise  with  10 
Percent  Uncorrelated  Noise  Added,  SNR  =  0.41 
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Figure  11-16.  Mean  and  Standard  Deviation  Vs  k,  Directional  Noise  with  10 
Percent  Uncorrelated  Noise  Added,  SNR  =  0.  83 
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Figure  11-17.  Mean  and  Standard  Deviation  Vs  SNR  at  k  =  0.  0  cycle/km,  Direc¬ 
tional  Noise  with  10-Percent  Uncorrelated  Noise  Added 
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Figure  11-18.  Mean  and  Standard  Deviation  Vs  SNR  at  k  =  0.  0  cycle/km,  Direc 
tional  Noise  with  20-Percent  Uncorrelated  Noise  Added 
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SECTION  III 

INTERPRETATION  OF  RESULTS 


In  all  work  discussed  in  this  report,  the  signal  was  considered 
to  have  infinite  velocity.  Figures  II- 1  through  II - 5  show  results  obtained  using 
the  13-element  array  conposed  of  the  inner  13  subarray  sites  of  the  Montana 
LASA  and  a  noise  model  consisting  of  completely  uncorrelated  noise. 

Data  presented  in  Figures  II- 1,  II-2,  andII-3  are  useful  in 
evaluating  the  relative  merits  of  the  conventional  and  high-resolution  schemes 
for  detection.  Here,  the  problem  is  to  select  a  threshold  level  which  will 
almost  always  be  higher  than  the  largest  peaks  resulting  from  the  ambient  noise. 
The  question  is  how  large  the  SNR  must  be  to  assure  that  the  signal  peak  will 
exceed  the  threshold. 

Note  that  arbitrary  scale  factors  were  applied  to  the  quantities 
plotted  in  theso  three  figures,  thereby  resulting  in  apparent  inconsistencies 
between  the  figures.  The  quantities  displayed  vary  by  orders  of  magnitude, 
and  scaling  was  necessary  to  <  i  oduce  meaningful  figures.  Since  the  important 
quantity  is  the  relationship  k.  tween  the  mean  and  the  standard  deviation  for 
a  given  processing  scheme,  the  same  scale  factor  was  applied  to  these  two 
measures. 

Figures  II- 1  and  II - 2  show  that  the  mean  and  standard  deviation 
for  both  schemes  have  local  maxima  near  0.  10  cycle/km;  the  spectral  window 
of  the  array  happens  to  have  a  peak  in  this  region.  Between  0.  02  and  0.  08 
cycle/km,  the  means  and  standard  deviations  are  quite  constant.  Moreover, 
the  levels  in  this  region  are  almost  exactly  the  same  as  those  observed  at  all 
wavenumbers  in  the  absence  of  signal.  Thus,  the  values  of  the  mean  and 
standard  deviation  at  0,  04  cycle/km  can  be  used  as  measures  of  the  behavior 
of  these  techniques  when  uncorrelated  noise  alone  is  present. 


Ill-  1 


science  services  division 


In  Figure  II-2,  the  mean  and  standard  deviation  for  the  con¬ 
ventional  spectrum  at  0.  04  cycle/km  are  0.  03  and  0.  Q25,  respectively.  If 
a  threshold  two  standard  deviations  above  the  mean  is  chosen,  its  value 
will  be  0.  08.  At  the  signal  location,  the  mean  and  standard  deviation  are 
0.  38  and  0.  14.  Thus,  the  mean  value  at  the  signal  location  is  slightly  more 
than  two  standard  deviations  above  the  threshold  level.  A  similar  calculation 
for  the  high-resolution  technique  shows  the  signal  mean  to  be  about  1.  8 
standard  deviations  above  the  threshold. 

In  Figure  II- 1  (SNR  =  0.  72),  both  techniques  result  in  signal 
means  that  are  about  1.  5  standard  deviations  above  the  threshold  (two  stand¬ 
ard  deviations  above  the  noise  mean  at  0.  04  cycle/km).  .  It  appears  that  the 
high- re  solution  scheme  offers  better  detectability  at  low  SNR's,  while  the 
converse  is  true  for  SNR's  above  about  0.  7. 

This  tendency  is  also  evinced  in  Figure  II-3,  where  the  abscissa 
measures  SNR.  Here,  the  standard  deviations  of  the  two  schemes  have  fairly 
similar  slopes,  but  the  slope  of  the  mean  for  the  conventional  spectrum  is 
appreciably  greater  than  that  for  the  high-resolution  spectrum.  This  is  con¬ 
sistent  with  the  superiority  of  tho  conventional  spectrum  at  higher  SNR's. 

In  summary,  for  an  SNR  of  1.0,  it  is  possible  to  choose  a 
threshold  level  which  is  two  standard  deviations  above  the  mean  level  for 
uncorrelated  ambient  noise  and  two  standard  deviations  below  the  mean  level 
for  the  signal  as  seen  by  the  conventional  spectrum.  While  slightly  inferior 
at  this  SNR,  the  high-resolution  scheme  apparently  becomes  somewhat  superior 
for  SNR  values  below  0.  7. 

Figure  H-4  illustrates  the  performance  of  these  processing 
schemes  for  locating  a  signal  in  uncorrelated  ambient  noise.  An  infinite - 
velocity  signal  was  embedded  in  each  of  100  independent  random-noise  vectors. 
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For  each  resultant  signal-plus-noise  vector,  the  conventional  spectrum  was 
evaluated  at  25  equispaced  locations  between  0.  0  and  0.  096  cycle/km  along 
the  N-S  axis,  and  the  peak-power  location  was  determined.  For  uncorrelated 
noise,  these  peak-power  locations  are  necessarily  the  same  as  those  resulting 
from  the  high-resolution  schemes  or  from  the  probabilistic  processor  (de¬ 
scribed  in  Appendix  A).  Figure  II-4  shows,  for  various  SNR,  the  number  of 
times  out  of  100  that  the  peak  power  occurs  at  each  of  these  wavenumber 
locations.  To  add  perspective,  note  that  at  1.  0  Hz  the  wavenumber  space 
distance  between  0.  044  and  0.  048  cycle/km  corresponds  to  about  700-km 
separation  between  epicenters,  while  the  distance  between  0.  068  and  0.  072 
cycle/km  corresponds  to  about  625-km  epicenter  separation.  Also  note 
that  the  locations  at  which  the  spectrum  is  evaluated  include  the  true  signal 
locations  but  only  a  small  fraction  of  the  total  number  of  other  locations  de¬ 
fined  by  a  rectangular  grid  having  a  0.  004- cycle/km  spacing  between  nearest 
locations.  Thus,  the  results  in  Figure  II-4  provide  an  upper  bound  on  the 
relative  fraction  of  correct  locations  that  can  be  expected. 

Subject  to  this  limitation,  it  is  seen  that  a  large  number  of 
incorrect  locations  occur  for  an  SNR  of  0.  25  or  0.  50.  At  an  SNR  of  1.  0,  how¬ 
ever,  the  signal  was  located  properly  98  times  out  of  100. 

When  a  velocity -preserving  stack  over  five  independent  fre¬ 
quencies  is  formed,  as  described  in  Section  H,  subsection  C,  significantly 
better  performance  results.  Figure  II - 5  shows  that  in  this  case,  the  0.25 
SNR  still  leads  to  a  number  of  incorrect  locations  but  for  an  SNR  of  0.  50  or 
1.  0,  the  signal  is  always  located  correctly.  The  data  in  Figure  II-5are 
appropriate  only  for  the  conventional  spectrum  technique.  Stacking  of  the 
other  types  of  spectra  would  lead  to  other  —  but  probably  similar  results, 
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implying  that  stacking  over  frequencies  in  the  frequency  domain  is  significantly 
superior  to  location  at  a  single  frequency.  When  stacking  is  done,  it  appears 
possible  to  properly  locate  the  signal  for  an  SNR  greater  than  0.  50.  It  is 
likely  that  the  analogous  procedure  of  selecting  the  largest  of  a  grid  of  time- 
domain  beamsteers  would  lead  to  similar  results. 

Figures  II-6,  II - 7,  and  II-8  illustrate  the  performance  of  a 
25-element  conventional  subarray  with  a  noise  model  that  is  largely  correlated. 
The  model  consists  of  an  isotropic  disk,  centered  at  the  origin  with  a  radius 
of  0.  1  cycle/km,  to  which  20-percent  uncorrelated  noise  has  been  added.  In 
this  case  and  all  subsequent  correlated  noise  cases,  the  procedure  consists 
of  generating  the  100  signal-plus-noise  vectors,  calculating  the  conventional 
and  high-resolution  spectra  at  11  equispaced  points  along  the  N-S  axis  for 
each  of  the  vectors,  normalizing  the  result  for  each  vector  so  that  the  11 
values  calculated  sums  to  unity,  and  then  calculating  the  mean  and  standard 
deviation  over  the  100  normalized  sets  of  data. 

Figure  II -6  shows  the  means  and  standard  deviations  for  con¬ 
ventional  and  high-resolution  spectra  as  a  function  of  wavenumber  at  an  SNR 
of  0.  0042.  This  small  amount  of  signal  has  virtually  no  effect  on  the  resultant 
spectra,  so  the  results  illustrate  the  appearance  of  noise  alone.  The  results 
for  an  SNR  of  0.  6  (Figure  U-7)  shows  that  it  is  not  possible  to  pick  a  thresh¬ 
old  for  either  spectral  technique  that  would  almost  always  be  above  the 
noise  peak  and  below  the  signal  peak.  Detection  with  a  standard  subarray  in 
this  type  of  highly  correlated  noise  would  be  possible  only  for  SNR's  much 
greater  than  0.  6.  The  data  in  Figures  II -6  and  II -7  indicate  that  the  high- 
resolution  technique  probably  can  detect  smaller  signals  than  can  the  con¬ 
ventional  technique.  This  conclusion  seems  consistent  with  the  data  of  Figure 
II-8,  where  the  means  and  standard  deviations  are  plotted  as  a  function  of 
SNR  at  the  signal  location,  0.  0  cycle/km. 
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The  noise  model  used  for  Figures  II - 9 »  JI-10,  and  II- II  consists 
of  a  disk  centered  at  0.  05  cycle/km  on  the  N-S  axis  with  0.  01-cycle/km 
radius  and  I-percent  uncorrelated  noise  added.  Figure  II -9  shows  the  means 
and  standard  deviations  as  a  function  of  wavenumber,  while  Figures  II- 10 
and  II- II  show  the  means  and  standard  deviations  plotted  vs  SNR  at  0.  0  and 
0.  05  cycle/km,  respectively.  For  ♦‘he  conventional  spectrum,  a  threshold 
at  0.  12  is  almost  always  above  the  largest  noise  peak.  Such  a  threshold, 
however,  is  above  the  mean  value  at  the  signal  location  for  any  of  the  SNR's 
considered.  An  appropriate  threshold  for  the  high-resolution  spectrum 
would  be  about  0.  095.  At  an  SNR  of  0.  4,  the  threshold  is  below  the  mean  at 
the  signal  location,  but  the  difference  is  less  than  one  standard  deviation.  The 
separation  increases  slightly  with  SNR  but  is  still  less  than  one  standard 
deviation  at  the  largest  SNR  considered,  0.  9.  The  conclusion  apparently  is 
that  for  subarray  processing  with  the  noise  field  consisting  of  very  concen¬ 
trated  bodywave  noise,  the  conventional  spectrum  has  questionable  value  for 
detection,  while  the  high-resolution  technique  has  marginal  value. 

Shown  in  Figures  11-12  through  II- 18  tre  means  and  standard 
deviations  resulting  from  a  noise  model  comprising  a  disk  with  a  0.  025- 
cycle/km  radius  centered  at  0.  05  cycle/km  on  the  N-S  axis.  Figures  11-12 
and  11-13  result  from  noise  models  having  1-percent  uncorrelated  noise 
added.  Figures  II- 14,  11-15,  11-16,  and  II- 17  have  10-percent  uncorrelated 
noise  added  to  the  noise  model.  Figure  11-18  has  20-percent  uncorrelated 
noise  added.  This  model  represents  bodywave  noise  which  is  mostly  direc¬ 
tional  but  fairly  widespread.  The  array  configuration  consists  of  the  center 
and  outer  15  locations  of  a  standard  LASA  subarray. 

Although  not  shown  in  these  figures,  the  data  for  all  three  per¬ 
centages  of  uncorrelated  noise  show  that  the  largest  peaks  —  when  noise 
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alone  is  present  —  occur  at  0.  0- cycle/km  wavenumber  for  both  spactral 
techniques.  Thus,  from  Figures  II- 13,  11-17,  andII-18,  the  means  and 
standard  deviations  at  the  origins  may  be  used  to  select  appropriate  thresh¬ 
old  levels.  If  the  threshold  is  chosen  to  be  two  standard  deviations  above 
the  noise  mean  for  the  conventional  spectrum,  the  resultant  level  is  above 
the  mean  at  the  signal  location  for  all  SNR's  considered.  For  the  high-res¬ 
olution  technique,  the  resultant  level  is  below  the  mean  at  the  signal 
location  for  SNR's  above  0.4;  but  in  all  cases,  the  difference  is  less  than  one 
standard  deviation.  The  conclusion  is  the  same  as  that  for  very  concentrated 
bodywave  noise:  the  conventional  spectrum  seemingly  has  little  detection 
potential,  while  the  high- resolution  technique  is  marginal  for  SNR's  above  0.4. 

Location  ability  using  the  probabilistic  processor  is  contrasted 
to  that  for  the  conventional  and  high-resolution  spectral  techniques  in  Figures 
11-19  and  11-20.  The  center  and  outer  15  locations  of  a  standard  LASA  sub¬ 
array  provided  the  16  elements  used,  and  10-percent  uncorrelated  noise  was 
added.  For  these  figures,  the  processor  outputs  were  calculated  at  13  wave- 
number  locations  with  a  0.  01-cycle/km  separation.  As  before,  the  results 
obtained  with  the  high- resolution  spectrum  are  identical  to  those  obtained  with 
the  conventional  spectrum.  The  figures  show  that  even  for  an  SNR  as  high  as 
6.  48,  any  of  the  spectral  techniques  applied  to  an  array  the  size  of  a  standard 
LASA  subarray  result  in  a  fairly  large  number  of  incorrect  locations. 
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APPENDIX  A 


DERIVATION  OF  PROCESSING  EQUATIONS 
A.  PROBABILISTIC  PROCESSING  EQUATION 


Frequency-domain  probabilistic  processing  is  a  method  using 
seismic  array  data  to  detect  plane-wave  signals  in  the  presence  of  ambient 
seismic  noise.  This  method  is  based  on  the  assumption  that  the  array  data 
is  Gaussian  with  zero  mean  and  known  crosspower  matrix  |^nJ  or  |q 
depending  on  the  absence  or  presence  of  signal.  The  array  output  is  con¬ 
sidered  to  be  k  complex  numbers  resulting  from  the  collateral  Fourier 
transform  of  time  gates  from  k  channels  of  simultaneous  array  time-domain 
data. 


The  criteria  for  estimating  the  presence  or  the  absence  of  a 
plane-wave  signal  is  analogous  to  testing  the  simple  hypothesis  that  the  ob¬ 
served  data  have  a  crosspower  matrix  JqJ ,  with  the  alternate  multiple 
hypotheses  that  the  data  are  from  a  set  of  crosspower  matrices  of  the  form 
[^s+n]’  suPerscr*Pt  m  re;fers  to  a  particular  plane-wave  signal  model. 

Derivation  of  the  Gaussian  multichannel  sampled  data  pro¬ 
cessing  equation  starts  with  the  multivariate  Gaussian  probability  distribution. 


P(xr  x2 . xk)  =  k  1-~  exp  J-)X|*T[x]  '  1  lx}  .  |  (A-l) 


In  this  equation,  matrix  notation  is  used  as  follows: 


)X| .  or  X. 
[x]  or  X 


=  i  set  of  complex  data  values  x  ,  x^ . x^» 

withXj  being  the  Fourier  transform  of  the  first 
time  gate 


*  crosspower  matrix  of  the  data;  i.  e.  , 


X..  =  x*x. 
lJ  1  J 
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The  superscript  T  designates  the  matrix  transpose;  the  superscript  -1  desig¬ 
nates  the  inverse  matrix;  and  the  superscript  *  signifies  the  matrix  complex 
conjugate.  Vertical  lines  bordering  a  square  matrix  indicate  the  matrix 
determinant. 


In  Equation  A-l,  k  is  the  number  of  data  values.  In  a  multi¬ 
channel  problem,  k  would  usually  be  equal  to  the  number  of  channels. 

The  detection  of  a  plane-wave  signal  in  the  presence  of  ambient 
seismic  noise  can  be  formulated  as  a  multipossibility  multichannel  Gaussian 
problem.  That  is,  one  can  assume  that  the  data  X^  is  a  sample  of  a  stationary 
Gaussian  time  series  have  a  k  by  k  crosspower  matrix  composed  of  noise 
[Q  1  ^r  one  of  a  set  of  signal-plus -noise  matrices  fcj 

To  completely  specify  the  problem  using  a  Bayesian  approach, 
the  a  priori  probabilities  P(Qn),  Possibilities  must  be  known 

or  estimated.  The  problem  is  to  obtain  the  a  posteriori  probabilities  of  noise 
alone  P(Qn  |X^)  and  of  signal  plus  noise  p(Q^.n  I  3C)  by  analyzing  the  data  X^. 

It  is  understood  that  Jq^J  ranges  over  all  allowable  plane-wave  signal 
models. 


In  the  seismic  detection  problem,  stands  for  ambient 
seismic  noise  and  P(Q  )  is  the  fraction  of  time  that  no  signal  is  present. 


m 


n 


Q‘*'  stands  for  a  particular  signal  model  in  the  presence  of  noise,  and 
s+n 

P(Q  )  is  the  fraction  of  time  that  this  signal  is  present, 
s+n 

A  statement  of  the  a  posteriori  probability  using  Bayes  theorem 
for  noise  alone  is 


p(Qn  |  X.)  = 


P(X:  |  Q_) 


n 


P(X.) 


P(Qn) 


(A-2) 
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and 


P(Qm  |  X.) 
s+n  1  i 


p(°r+a> 

nki°r+j 


-<T  [C]  xii 


(A-9) 


P(X.) 


The  denominator  of  Equations  A -8  and  A-9  normalizes  P(Q  X.) 

■*  n  i 

and  P(Qm  lx.)  so  that  their  sum  is  1.  Equation  A-9  may  be  normalized  by 
s+n  '  i 

replacing  the  denominator  in  Equation  A-9  by  the  summation  of  the  numerator 
over  the  possible  values  of  m.  In  this  case,  the  possibility  that  no  signal 
is  present  is  removed.  Furthermore,  if  it  is  assumed  that  the  probabilities 
of  the  occurrence  of  the  various  signals  is  equally  likely,  Equation  A-9  becomes 


m  I<CJ 

P(Q“  I  x.)  =  — 


exp  J-X 


^lQm  I 

m  s+n1 


r  fa]"*, | 

“p  KT  fcj ' x. 


(A- 10) 


For  locating  a  plane -wave  in  noise,  the  .crosspower  matrix 


becomes 


[<|  =  [n + vv*T] 


where  N  is  the  noise  crosspower  matrix  and 


ik  •  r. 


ik  •  r 


V  = 


ik  •  r 
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with 


k  =  vector  wavenumber 
th 

r  =  p  array  coordinate  corresponding  to  the 
^  p  data  channel 

Thus,  the  superscript  m  denotes  a  plane-wave  signal  located 
at  one  of  the  possible  points  in  K-space. 


matrix  may  be  written  as 


For  computational  purposes  the  inverse  of  the  crosspower 
* 


£n  +  vv*T] 


-1 


=  N 


- 1  *T  - 1 
-1  N  VV  N 

*T  - 1 

1  +  V  N  V 


(A  - 1 1 1 


B.  HIGH-RESOLUTION  FREQUENC Y-WAVENUMBER  SPECTRA 


High-resolution  frequency-wavenumber  spectra  are  computed 
from  a  frequency-domain  multichannel  filter.  In  the  filter  synthesis,  the 
signal  crosspower  matrix  is  assumed  to  be  uncorrelated  Gaussian  noise,  and 
the  noise  crosspower  matrix  is  measured  data. 

The  multichannel  filter  designed  at  a  single  frequency  is 


WI  +  x.x  T 

1  1 


(A-12) 


— - -  "  *pr  ' 1  *T  "  *T 

*The  formula  used  is  (A  +  XY  )  =  R-RX(1  +  X  RY)  Y  R,  where 

- 1  *T. 

X  and  Y  are  row  matrices  (vectors)  and  R  =  A  .  Both  A  and  (A  +  XY  ) 

are  assumed  to  be  nonsingular. 
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where 


fol 


* 

F  = 


n 


> 


r  = 

n 


> 


V  J 


l°J 


k  and  X.  are  defined  in  the  previous  sections;  w  is  a  scale  factor;  and  I 
is  the  identity  matrix.  The  subscript  n  (1  £  n  £k)  specifies  the  reference 
channel  used  in  the  filter  design. 

th 

Using  the  n  channel  as  reference,  the  solution  to  Equation 

A- 12  is 


F 


* 


n 


1 

(w  +  X*TX.) 

i  i 


r 

n 


(A  -  13) 


A  high-resolution  fr equency-wa venumber  spectra  may  be 
computed  from  Equation  A-13  with  a  single  sensor  as  reference,  or 


*  T 
F  F 
n  n 


(A -14) 


Summation  of  the  responses  from  using  each  channel  as  re 
ference  results  in 


T  *  T 
V  F  F 
n  n 


u  n 


(A  - 15) 
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Letting  B  =  . 

&  n  N  n 


Equation  A-  14  becomes 


T  *  T  * 
V  F  F  V 

n  n 


(A  - 16) 


±rT 

w  +  X.  T  X. 

i  i 


* 


(w  +  X*TX.) 


VTX.X*T[yB*BTl 

i  i  /  ^  n  n 

L  n  ■* 


However,  when  all  sensors  are  used  as  reference, 


z 


*  T 

B  B  =  I 
n  n 


n 


and  Equation  A- 16  simplifies  to 


T  *  x  *  T  *  *  (2w  +  x  X  ) 
V  V  -  (V  X.  )(V  X.  )  - - - — 

'  (w  +  X*TX.) 

I  1 


(A- 17) 
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APPENDIX  B 

DETERMINANT  EVALUATION 


It  is  desired  to  develop  an  efficient  method  for  the  evaluation 

•  *T  i  -* 

of  I  N  +  VV  I  for  a  large  number  of  values  of  k.  Here,  N  is  a  covariance 

matrix  and  V  =  V(k)  is  a  unit  plane -wave  transform  for  a  signal  at  the  vector 

~4 

k  in  wavenumber  space. 

First,  the  following  theorem  is  proven.  Let  I  denote  the 
T  n 

(n  x  n)  identity  matrix  and  W  =  (w , ,  w_ ,  .  .  .  w  )  denote  an  n-dimensional 

n  1  2  n 

complex  vector.  Then, 

1 1  +  W  W*T  1=1  +  w*Tw 

1  n  n  n  1  n  n 

This  is  proved  by  induction  on  the  integer  n. 

Clearly  the  theorem  holds  for  n  =  1.  Assume  that  it  holds  for 
n  =  r.  The  problem  is  to  show  that  the  theorem  holds  for  n  s  r  +  1.  Let 


Z  |  =  [i  +  w  w*T" 

r  '  r  r  r 


Partition  Z  ,  such  that 
r+1 


M  ■ 


1  +  w  .  .  w  .  . 
r+1  r+1 


where 


MT  =  (wlwr+l . wrWr+l) 


Leta  =  1  +  w  ,w  ,•  Then  I Z  ,  I  can  be  written  as 
r+1  r+1  1  r+J  ' 
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1  *T  i 

|zr+1l  ■  lzt-M“  M  M« 


[zr  -  =|ir  +  a'1wrw*Tj 


which  may  be  verified  by  direct  computation.  Then  by  assumption, 


Thus, 


I  zr+i  I  ~  11,  +  a'1  WrW*T|.|a| 


_  l  jjj'p 

a  (]  +  a  w  W  )a 

r  r 


I  Zr+1  I  " 


♦T 

a  +  W  W 
r  r 


a  i  +  W  -  W  ‘  +  W  TW 

r+ 1  r+ 1  r  r 

a  1  +  W**  W 

r+ 1  r+ 1 


*T  r 

To  evaluate  N  +  VV  |,  consider 


Q.E.D. 


N  +  VV 

Since  N  is  positive  definite  Hermitian,  there  exists  a  matrix  S  such  that 

*T 

SNS  =  I 


It  follows  that 


*T  *T  *T  *T  *T  *T 

S(N+VV  )  S  a  SNS  +  SVV  S  a  I  +  WW 
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where  W  =  SV.  Then 

|  S(N  +  VV*T)  S*T  !  =  |  S  |  I  N  +  VV*T  |  |  S*T  |  =  1 1  +  ww*T  | 

Making  use  of  the  previously  proven  theorem,  this  becomes 

3#crr 

In  +  vv  |  *  |  n  |  (i  +  w  w) 

i  *T  i 

which  is  an  efficient  method  for  evaluating  |  N  +  VV  | . 
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Various  noise  models,  made  up  of  mostly  organized  noise,  were  used.C^our  spectral  tech¬ 
niques  —  the  conventional  spectrum,  two  high- re  solution  spectra,  and  probabilistic  pro¬ 
cessor —were  evaluated.  When  wavenumber  spectra  at  a  single  frequency  are  used  to  locate 
a  signal  In  uncorrelated  noise,  all  four  techniques  must  yield  exactly  the  same  location.  It 
was  found  that  events  with  a  SNR  of  1. 0  or  greater  were  almost  always  located  properly  within 
the  resolving  power  of  the  grid  points  at  which  the  spectra  were  calculated.  Wideband 
frequency-domain  location  can  be  accomplished  by  computing  wavenumber  spectra  at  several 
independent  frequencies  and  summing  the  spectra  in  a  velocity-preserving  manner  before 
choosing  the  peak-power  location.  Conventional  spectra  of  signals  in  uncorrelated  noise  were 
computed  at  five  independent  frequencies  and  summed.  This  procedure  results  in  no  incorrect 
locations  for  SNR'  s  of  0.  5  or  greater,  thereby  indicating  its  significant  superiority  over  the 
location  ability  of  single -frequency  spectra.  The  wideband  wavenumber  spectral  technique  is 
roughly  analogous  to  time-domain  beamforming,  and  the  results  obtained  here  may  be  indica¬ 
tive  of  those  achievable  in  the  time  domain.  The  performance  of  the  small  arrays  in  organized 
noise  was  found  to  be  markedly  inferior  to  the  performance  of  the  large  array  using  uncorre¬ 
lated  noise.  The  conventional  spectrum  seems  to  have  no  practical  detection  ability,  while  the 
high-resolution  techniques  are  perhapB  marginal  in  this  respect.  All  techniques  located  a 
significant  number  of  events  improperly,  even  when  the  SNR  was  as  high  as  6.  48. 
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